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It is widely known that the paradigmatic Chirikov-Taylor model presents enhanced diffusion for specific 
intervals of its stochasticity parameter due to islands of stability, which are elliptic orbits surrounding accele- 
rator mode fixed points. In contrast to the normal diffusion, its effect has never been analytically calculated. 
Here we introduce a differential form for the Perron-Frobenius evolution operator in which normal diffusion 
and superdiffusion are treated separately through phases formed by angular wavenumbers. The superdiffusion 
coefficient is then calculated analitically resulting in a Schloemilch series with an exponent /? = 3/2 for the 
divergences. Numerical simulations support our results. 

PACS numbers: 05.45.-a, 05.20.-y, ()5.60.-k 



One of the most fundamental questions in the study of 
transport theory is to understand strange diffusion, whereby a 
mean-squared displacement follows an asymptotic time power 
law (a;^) ~ t'^ for/3 ^ 1, starting from nonlinear microscopic 
equations of dynamics. The well-known Brownian (normal) 
diffusion is characterized by the exponent /3 = 1, while the 
strange diffusion is subdivided into two types: superdiffusive 
(/3 > 1) and subdiffusive (/3 < 1). The strange transport has 
been studied and observed in a wide range of systems such 
as plasmas fl], flows and turbulence [2], amorphous semi- 
conductors, porous media, glasses, granular matter [3], one- 
dimensional intermittent maps ||4|], and biological cell motility 

fl. 

In this Letter, superdiffusion is studied analytically in the 
paradigmatic Chirikov-Taylor standard map: 



In+l 
0n+l 



In + K sin Or, 



I, 



n+l 



mod 27r, 



(1) 



defined on the cylinder —Tr<0<'K,~oo<I<oo, where 
K is the stochasticity parameter In the past, many authors 
investigated the diffusive behavior of chaotic systems on the 
basis of map ([T]i through the diffusion coefficient D, defined 



D = 



{{in-ior) 



./3-1 



2n 



n> 1. 



(2) 



Brownian diffusion emerges when the last KAM barrier is de- 
stroyed at the critical parameter value Kc ~ 0.9716 la]. The 
corresponding diffusion coefficient D{K) exhibits an oscil- 
latory feature for K > K^. This behavior was first numer- 
ically observed by Chirikov \m and subsequently calculated 
by Rechester et at. [8'|. In addition to the normal diffu- 
sion, the map ([T]l also exhibits superdiffusion of the Levy type 
(1 < /? < 2) due to accelerator mode islands, which are cre- 
ated around a stable Q-periodic orbit {(!*, ^r*)}i<n<0' satis- 
fying LSJ 



i:.+Q-i:^2ln, i^£sin0; 



2/7r, 



(3) 



for any integer I ^ 0. For Q 
stability condition jlOll : 



1 we have the following 



Ss : |2 + /\cos6'| < 2. 



(4) 



Typically, trajectories diffuse normally, although some of 
them may be dragged along the accelerated modes. Such co- 
existence between chaotic and regular orbits results in a com- 
plex dynamics whose exponent is certainly less than the bal- 
listic value P — 2. This fact has motivated the use of suit- 
able statistical theories such as continuous-time random walk 
(CTRW) LUJ and fractional kinetics (FK) [12] in the numer- 
ical study of the superdiffusion. Here, the superdiffusion co- 
efficient and its corresponding transport exponent are analyti- 
cally calculated from first principles, i.e., from dynamics ruled 
by Eq. dTJ. 

The analysis of the map ^ is best carried out in Fourier 
space. The Fourier expansion of conditional probability den- 
sity pn, in which an initial state (Iq, 9q) at n = evolves to a 
final state (/, 6) at step n, can be written in the form 



,(/,0|/o,0o)-^e™'an(m,/). 



(5) 



Relations between normal diffusion and Pollicott-Ruelle res- 
onances for a general class of area-preserving maps that in- 
cludes Eq. ([T]i as a particular case was studied in Ref. OlSll 
through the evolution of transformed amplitudes an{m, q) = 
(27r)"i/d/a„(m, /)e"*'?'^. On the other hand, as it will 
be seen latter, the ampltudes an{iTi,I) are proving to be 
more suitable in the description of superdiffusion. The dis- 
crete time evolution of these amplitudes are also ruled by the 
PeiTon-Frobenius (PF) operator U, defined as dn{iTi,I) = 
C/a„_i(m, /). More explicitly, starting from the evolution of 
an{m, q) we have 



an{m,q) == (27r)" 



d9dl e 



-i(me+ql) 



n=l 



X I de'di'5{i -I' -K sin e')5{e -9' -I) 

xY^Jdq' £'('"''''+'' ^') a„_i(m', q'). (6) 



The useful PF operator for a„(m, /) can be obtained in three 
steps: i) integrating Eq. ^ on 9 and /'; ii) multiplying the 
two sides of the resulting equation by exp(jgJ') and integrat- 
ing both on q; and Hi) considering the translation operator 
a{I + x) = ex-p{xdi)d{I). Finally, the dynamics of the sys- 
tem ([T]| induces the following law of evolution 



a„(m,/) = e-'™^^J„ 



,{iKdi)dr,-i{m',I), (7) 



where Jm{x) is the Bessel function of the first kind. Here- 
after, we will use the abbreviation Jm{iKdi) = Jm- 
The initial probability amplitude is given by ao{m,I) = 
(27r)^^e^""'^°(5(/ — /q), which represents the certainty of an 
arbitrary initial condition I — Iq and 9 — 9o at n ~ 0. 
Besides, we must define the density to generate initial val- 
ues (/q, 9o), which are independent from the dynamics. Usu- 
ally, the initial values are assumed as uniformly distributed in 
their respective domains jl41]. In this sense, numerical calcu- 
lations of (|2]) are performed as Dnum — N^^ J2i=i{-^n,j — 
Iq,jY /2n, where iV ^ 1 is the number of random initial con- 
ditions generated for each value of K. 

In order to calculate the superdiffusion coefficient Dsup it 
is useful to calculate the normal diffusion coeffcient Dnor by 
means of this new approach. In this context, arbitrary mo- 
ments depend only on the state m = at instant n: (Jp) = 
{2Tr)^[{idg)PaniO,q)]g=o = (2^) / d/5„(0, 7)/^. Thus, we 
can decompose the projected PF operator PC/" into two parts: 
PU" = PU"P + PU"Q, where P and Q ai-e mutually or- 
thogonal projection operators that represent the states m — Q 
and m ^ Q, respectively. On the other hand, PU^Q can be 



neglected due to ao(m ^ 0,q) 



oc e 



-imOo 



whose expected 



value disappears at random initial conditions on [— tt, tt). 



Hence, PU^^P rules the diffusion in the normal regime 1113 1. 
Now consider the leading case where mi = ra2 — . . . ~ 
nin = 0, which gives (l" = Jp" = 1 + n{K'^/^)d'j + 0{dj) 
and the quasilinear diffusion Dqi = K'^/4. Corrections to 
this value are obtained by the introduction of a small fc-piece 
of wavenumbers in a quasilinear stripe of length n, as shown 
in the example below for fc = 3: 



Thus, in order to obtain only relevant contributions to the 
diffusion coefficient, we should perform summation in Eqs. 
( 19110b over all integers subject to the constraint nii + m2 + 
7713 = 0. Besides, Bessel operators follow the same series ex- 
pansions of Bessel functions, namely J„i oc 9™[1 + 0{dj)]. 
Evidently, wavenumbers that contribute with terms 0{dj) can 
be neglected in the calculation of the moment (/^) resulting, 
in addition to ms = —{mi + 7712), in the three independent 
j-constraints |77ii| + \^j77ii + aj77i2\ + \77i1 + 7712! < 2, 
where (/ij-jCTj) — {(0, 0), (1, — 1), (1, 2)}. Applying these 
constraints, the sum V^ ^ ^ f^ becomes Vti + CI2 where 



hi = 2J2{K)JoJl 

h2 = 2Jl{K)Jl + AJ2{2K)Ji{2K)jj2 + .h{K) 
X [-h{2K) + J:,{K) + Ji{2K) - Ji{K)]JI 

Finally, the normal diffusion coefficient is calculated as 



(11) 



(12) 



^"""^ = ^ I dI[l+77{(ll+VL2)]{D-i^+77dl)6{I)l' 



Dqi 2n , 

= 1 - 2J2(^) - 2Ji (X) + 2Jl{K) 



(13) 



for 77, 3> 1, where T = I — Iq. The terms of the series ex- 
pansion ( fljl l coincides with the Rechester et al. results ||8|], 
noting that -J2{2K)Ji{2K) w J^{K) for if > 1. Further 
corrections can be obtained by increasing the fc-piece. 

We turn now to the analysis of superdiffusion. In con- 
trast to the dynamics of normal process, accelerator modes 
allow the PF operator t/" to load phase terms of the 
type exp(7/^^ TTifc) without the need of previous condition 
Tlik "^fe ~ ^- ^^ ^^^ other hand, accelerator modes also con- 
tribute to the corrections of the quasilinear value, although 
much stronger than the quasilinear one. Therefore, the stripes 
of wavenumbers that represent the leading contributions for 
the superdiffusion are 

(0, 0) . . . (0, 0) (0, 777) (777, 777) .. . (777, 777), (14) 

71— r— 1 1 r 



(0, 0) . . . (0, 7ni)(777l, 7772)(7772, 7773)(7773, 0) . . . (0, 0) . (8) 
fc=3 

Considering k ^ nwe can neglect border effects, so that we 
have 77 configurations of the type (O that gives the same result 
for PU"P. Adding to these the quasilinear configuration we 
have 

an(0,/)-^(l + 77^^^f7j JS6iI-Io\ (9) 

\ rni m2 rris / 

where the operator fl is given by 

u L ' ry — jyy T '^ ^ 

X [Jmi-m2{'m2K) + Jmi-m2i'm3K) + Jmi-mJ 

X [Jm2-7n3{m3K) + J„i2-m3]'X Jma- (10) 



(rn,m) . . . {m,m), 



(15) 



corresponding to particular configurations of PU^Q and 
QU"'Q, respectively. Adding the contributions due to stripes 
( fT4l ) and ( fTSl ), the density p„ assumes the form 



1 " 
.iI,9\Io,9o) = ^EE 



Jm[{9-eo)-rI] 



r—1 m^O 

x<i>rMK)J'o^iI - lo) (16) 

where, neglecting 0{J^) terms and considering n ^ 1, we 
have 



J^rnimK), 

<i>rMK) - { 0{J^), 

J2'LiJoiJmK), 



2 < r < 77, (17) 



The leading terms in the Eq. ( fTSI l are phases restricted to the 
following conditions 



/ ^ /„ 



2pvr, 



(18) 



for any integer p, including /q = 2p7r as an alternative so- 
lution due to delta function 6{I — Iq). In other words, only 
initial conditions that start from the vicinity of /q = 2pTT (or 
final conditions that reach it) for any 6*0 G Se and finish as a 
Q-periodic orbit (1 <Q< n) must contribute to the superdif- 
fusive transport. This scenario includes not only accelerator 
modes, but also irregular trajectories that are dragged along 
them. Such combination of regular and irregular paths makes 
it difficult to know precisely the probability measure for an- 
gular variables {6, da). In this sense, the angular amplitude 
along the accelerated dynamics can be estimated as 



dedOo e""(^-^«) 



-^A^A^o = ^§^, (19) 



where the maximal width of stability A6'o is given by Sg'- 

A6lo = arccos(-4/i^) - tt/2 w i/K, (20) 

valid for K > 4. The averaged phase ( fT9] l is justified con- 
sidering that: i) its value is lower than 2A9A9o, where the 
factor 2 is due to the existence of two gaps; and ii) the same 
integral must be proportional to m^^. In Fig. [T]we illustrate 
this scenario for K = 1.1 x 27r (Q = 3 accelerated mode H), 
showing two gaps with widths A^o as foreseen by Eqs. (|4]l 
and (|20] |. Note that orbits lying in the gaps are captured by 
condition (fTsT l while the others develop Brownian motion. 

The stability interval for accelerated modes in terms of the 
parameter K is given by 



>K 



2\1\tt <K < 2\l\TTy/l + {2/lTTy 



(21) 



obtained through Eqs. (|3]l and (|4]l for Q = 1. The stability 
windows Sg and Sk are maximal since they contain all other 
Q > 1 corresponding stability intervals for each correspond- 
ing value of I IlZIIioI]. We can now compare the weight coeffi- 
cients ^i^jn{K) and $„ „i(if) subject to Sk- For K w 2\1\tt, 
the sum X^ms^o "^~^^i,m(^) can be estimated as follows; 



E 



J^rnimK) 



vir 



V m—1 



E(-ir^'c™-^F ^22) 



(m - l/2)-5/2 



m 



-5/2 



, noting that this se- 



where c„ 

ries is convergent by the alternating series criterium. On the 
other hand, ^n^m{K) can be assessed for n 3> 1 using a 
Schloemilch series [11511 ; 






j=i 



valid on the interval 



(23) 



21't: < \m\K < 2{l' + l)Tr, 
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FIG. 1: Measurements of 9 (n=100) versus 9o from numerical sim- 
ulations performed for iV = 5000, K = 1.1 x 27r, Jo = 27r, 
and 6o uniformly distributed on [—n,n). Accelerated orbits clus- 
ter along the diagonal 6 = 9q on the gaps located inside the interval 

7r/2 < |6lo| < arccos(-4/J-i') f» 2.188. 



where I' is a positive integer Comparing the restriction 
( |24| ) with the stability condition ( ITTT i it is easy to see that, 
in the limit ti ^ oo, the superdiffusion rate will diverge 
on i^ = 2\l\Tr by setting I' = |7tt,^|. Furthermore, due 
to the strong decay of angular amplitudes ( fT9] l, the lead- 
ing wavenumbers are m = ±1. Finally, the mean-squared 
displacement due to accelerated dynamics is calculated by 
^X'^(n)) = J dlpacc (2^, '^)2^'^ resulting for the diffusion 



' sup 



D, 



Dn 



16 

if 

16 

'k 



oo / n \ 

m=l y i=l / 



Y^MjK) for KeSK 



(25) 



3 = 1 



(24) 



and zero otherwise. 

We have made extensive numerical studies of the superdif- 
fusion coefficient ( |25] l, including the normal diffusion term 
( fT3] l. A particular case for n = 100 is showed in Fig. |2] 
It is clear that there is good agreement between theoretical 
and numerical results, especially when it is observed the ac- 
cumulation of enhanced occurrences to the right of theoretical 
spikes in contrast to the their quasi absence on the left. In ad- 
dition to the leading peaks located on the intervals pT|) with 
divergences at K ~ 2\1\7t for n -^ oo, map ([T]l has a more 
intrincated chain of secondary peaks whose phenomenology 
has been well studied in Ref. JmllTlllSll . In order to describe 



them, we need to retain more terms in the evaluation of Eq. 
(fTSI l including quasilinear stripes with a growing number of 
"failures" (sequence of different values for wavenumbers). 
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FIG. 2: Theoretical diffusion rate D/Dqi (solid line) compared with 
normal rate Dnor/ Dqi (dotted line) and measurements of D / Dqi 
from numerical simulations (A) for n — 100 and N = 5000. First 
and second theoretical diffusion peaks occur for values 14.63 (K — 
6.30) and 6.07 {K = 12.58), respectively. 

An important question concerns the determination of the 
superdiffusion exponent. As pointed out by Zaslavsky and 
Edelman [17], typically it is hardly possible to get a theoret- 



ical value of p. If the power law D 



sup 



i/3-l 



IS asymptot- 



ically valid, then /3 > 1 may exist only for values of K for 
which Dsup diverges. Using the asymptotic form of Bessel 
functions and noting that X]?=i J^^^^ ^ 2n^/^ for n 3> 1, 
we can evaluate Eq. dZSl ) for K = 2|Z|7r as 



D 



sup 



D„ 



|/|3/2 



(26) 



resulting (3 = 3/2. In Ref. llSll there is a comparison between 
CTRW and FK formalisms applied to the study of superdiffu- 
sion of the map (HJ. In both studies the numerical values of 
(3 for K ^ 1.03084 x 27: (Q — 5 accelerated mode) are also 
very close to 3/2, namely /3 = 1.42±0.15. Note that Eq. (|26] | 
not only represents a class of universality for /3 but also sug- 
gests the existence of other classes related to secondary diffu- 
sion peaks. In a similar way, further corrections for Dgup may 
give secondary divergence terms for tighter iiT-intervals such 
that, together, they form a multifractal layer for the mean- 



squared displacement in the form <J/^) 
where /3j — j3j (K) as numerically observed in 



E^/3. 



7,/5^. 
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